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We study the family of static and moving dark solitons in quasi-one-dimensional dipolar Bose- 
Einstein condensates, exploring their modified form and interactions. The density dip of the soliton 
acts as a giant anti-dipole which adds a non-local contribution to the conventional local soliton- 
soliton interaction. We map out the stability diagram as a function of the strength and polarization 
direction of the atomic dipoles, identifying both roton and phonon instabilities. Away from these 
instabilities, the solitons collide elastically. Varying the polarization direction relative to the con¬ 
densate axis enables tuning of this non-local interaction between repulsive and attractive; the latter 
case supports unusual dark soliton bound states. Remarkably, these bound states are themselves 
shown to behave like solitons, emerging unscathed from collisions with each other. 

PACS numbers: 03.75.-b,03.75.Lm,05.45.Yv 


I. INTRODUCTION 

Solitons - waves which propagate without dispersion 
- are a paradigm of nonlinear physics, occurring across 
systems as varied as water, optical fibres, spin chains, 
the human circulatory system and atomic Bose-Einstein 
condensates (BECs) [T]. Much of the soliton behaviour 
across these diverse systems, including their interactions 
and collisions, are universal [2]. The extreme controlla¬ 
bility of the nonlinearity, dimensionality and external po¬ 
tential in atomic BECs makes them an ideal playground 
for studying solitons Em, with experimental demon¬ 
strations of bright [MU, dark [EHIT] and dark-bright 
[muH] solitons. Moreover, these “matter wave solitons” 
provide insight into the interplay of solitons with quan¬ 
tum coherence [19], and have potential applications in 
interferometry [ZliniEolEi], surface force detection [22] 
and as robust quantum-information carriers [23] . 

In conventional media with local nonlinearity, soliton 
behaviour is well-established EIMI, e.g., solitons interact 
only at short distances (when their field profiles overlap). 
A current direction in soliton research is the addition 
of non-local nonlinearity. This enriches the soliton be¬ 
haviour, e.g., promoting non-local interactions between 
solitons, and offers prospects for studying wave analogs 
of particles with long-range interactions and emulating 
complex nonlinear networks [25]. To date, such “non¬ 
local solitons” have been observed in liquid crystals [26]- 
[28] . thermo-nonlinear optical fibres [25] and liquids [29] . 
and optical ring fibres [30]. In these cases, the non-local 
response is provided by molecular re-orientations, heat 
conduction and acoustic waves, respectively. Typically 
these nonlinearities lead to retarded interactions. 

The advent of BECs of atoms with sizeable mag¬ 
netic dipole moments - ^^Cr [31] [32], ^^^Dy [33l|34] and 
^^^Er [35] - allows for the study of superfluids with non¬ 
local nonlinearities [36], arising from the dipole-dipole 
(DD) atomic interactions. Remarkably, the ratio of local 
to non-local interactions can be directly tuned through 


Eeshbach resonances This has opened the door 

to experimental observations of magnetostriction [38] . 
anisotropic collapse [39] and self-organized droplet phases 
[34] . A host of predictions such as 2D bright solitons m 
and stabilised 3D dark solitary waves im await experi¬ 
mental verification. 

In this paper we obtain numerically the family of static 
and moving dark solitons in a quasi-ID dipolar BEG as 
stationary solutions in static and moving frames. We 
map out the stablity diagram and form of the soliton 
solutions as a function of the strength and polarization 
direction of the dipoles, covering experimentally-relevant 
parameters. The dark solitons can acquire dramatically 
modified profiles, including ripples close to the roton in¬ 
stability. Moreoever we show analytically that the ab¬ 
sence of atoms in a dark soliton causes it to act like a 
giant anti-dipole [42]. This induces a non-local contri¬ 
bution to the soliton-soliton interactions which scales as 
1/z^. This non-local contribution is effectively instanta¬ 
neous and experimentally tunable in both strength and 
sign. When the non-local interaction between the solitons 
is attractive, the balance with the conventional repulsive 
soliton-soliton interaction supports bound states of two 
dark solitons, for which no analog exists in conventional 
condensates. We show that these bound states them¬ 
selves have soliton-like properties in that they emerge 
unscathed from collisions with each other. 


II. THEORETICAL FRAMEWORK 

We consider a dilute, ultracold BEG of atoms with 
mass m and polarized magnetic dipoles. At such low 
energies the interactions between the atoms can be de¬ 
scribed using a universal pseudo-potential [36] . 

U{r-r')=gS{r-r')+UM{r-r'). (1) 

The first term is a contact interaction accounting for 
the van der Waals interactions characterized by coefii- 


2 


cient g. This gives rise to a local nonlinearity/mean- 
field potential ^n(r), where n(r) is the atomic density. 
The second term gives the dipole-dipole interactions, 
with t/dd(r ~ r') = Cdd(l ~ 3cos^ 0)/Att\y — r'p, where 
0 is the angle between the polarization direction and 
the inter-atom vector r — r', and Cdd (which is con¬ 
ventionally positive) characterizes the strength of the 
dipoles. This contribution to the interactions gives rise 
to a non-local nonlinearity/dipolar mean-field potential 
^(r) = /[/dd(i* — r')n(r') dr'. At the magic angle 
~ the dipole-dipole interactions disappear. For 
0 < Om the dipoles he dominantly head-to-tail and at¬ 
tract, while for 0 > 0^ they lie dominantly side-by-side 
and repel. However, by fast rotation of the polariza¬ 
tion direction the time-averaged dipole-dipole interaction 
can be effectively reversed (Cdd < 0)? in which case the 
dipoles repel when head-to-tail and attract when side-by- 
side [43] . 

In three-dimensional geometries, dark solitons are 
prone to transverse excitation of the nodal line, the 
so-called “snake instability”. However, in quasi-one- 
dimensional geometries this instability is prevented and 
dark solitons become long-lived 0 la m. We con¬ 
sider such a quasi-ID waveguide geometry, aligned along 
z. The confinement in the transverse directions is as¬ 
sumed to be harmonic with the form, l/(r) = mco’^r^/2, 
where u;± is the corresponding trap frequency. The quasi- 
ID limit is reached when this confinement is sufficiently 
strong {huj±_ where g is the BEC chemical potential) 
that the condensate wavefunction approaches the har¬ 
monic oscillator state in the transverse direction [45j [46] , 
with lengthscale l± = yjTijmu^. In this regime the 3D 
condensate can be parameterized via a ID mean-field 
wavefunction which obeys an effective ID Gross- 

Pitaevskii equation (GPE) [471148] . 

Here 4 >id {z^ t) = f Uid(z — z')\i^\^ dz' is the effective ID 
dipolar potential, with associated pseudo-potential [49] . 


Um{u)=Uo 


2rt—\/^(l+rt^)e^ /^erfc 




(3) 

where u = \z\/l^ and Uq = Cdd(l + 3cos26>)/327r/^. 
While the above mean-field model ignores finite temper¬ 
ature effects, it is worth noting that the salient physi¬ 
cal behaviour of the (non-dipolar) dark solitons observed 
experimentally to date are well-described by the zero- 
temperature mean-field Gross-Pitaevskii model 


III. STABILITY OF THE HOMOGENEOUS 
SYSTEM 

The stability of the homogeneous ID dipolar BEC has 
been established in Refs. gaiiH] for dipoles aligned along 



FIG. 1. (Color online) Stability diagrams in the 
(^, £dd)—plane for the homogeneous quasi-lD dipolar BEC 
for (a) ^ > 0 and (b) ^ < 0. Depicted are the regions of sta¬ 
bility (white), phonon instability (blue) and roton instability 
(red). The magic angle 0 ^ is highlighted. 


z {0 = 0); here we map out the full parameter space 
(0 < 6> < tt). The ground state has uniform ID density 
no- Its chemical potential (the eigenvalue associated with 
the right-hand side of Eq. Q) is /io = ^o^/27r/^ + 4>o, 
where the first term represents the van der Waals mean- 
field potential and 4>o = —(7dd^o[l + 3 cos 2^]/247r/^ rep¬ 
resents the dipole-dipole mean-field potential. The char¬ 
acteristic length and speed scales are the healin g length 
^ = h/y/mfio and speed of sound c = golm; a time- 
scale follows as r = ^/c. We parameterise the transverse 
BEC size via a = l±/^ where the quasi-lD approximation 
requires a ^ 1 [45]. We take the arbitrary value a = 0.2 
throughout this paper, but our qualitative findings are 
independent of a. 

It is convenient to parameterize the dipole strength via 
the ratio Sdd = Cdd/^g [36]. In the dipolar condensates 
produced to date, this ratio has the natural value Sdd = 
0.16 for ^^Cr [64], Sdd = 0-4 for ^^^Er [35] and Sdd = 1-4 
for ^^^Dy[33j|65]. However, g and Cdd can both be tuned 
in size and sign [36l [43] , and so we will typically consider 
a range of Cdd with both negative and positive g and 
Cdd (later, when considering soliton collisions and bound 
states, we will focus on the parameters for ^^^Er, and 
comment on the wider dependence on 5dd)- 

The homogeneous condensate is unstable for certain 
Cdd, g and e mm ; the stability in the (^, 5dd)-plane 
is depicted in Eig. [^for (a) ^ > 0 and (b) < 0. The 

condensate suffers two key instabilities - the phonon in¬ 
stability and the roton instability. The phonon instabil¬ 
ity refers to unstable growth of low /c-modes. It arises 
when the net interactions become attractive, i.e. when 
go < 0 (blue shaded regions). Consider, for example, 
conventional dipoles (Cdd >0): for ^ = 0, the phonon 
instability arises when the attraction of the end-to-end 
dipoles dominates the van der Waals repulsion, while for 
^ = 7r/2 it arises when the repulsive side-by-side dipoles 
are outweighed by attractive van der Waals interactions. 
The reverse is true for Cdd <0. In both cases a long 
wavelength collapse is induced. 

The dispersion relation can also feature a roton-like 
dip at finite momenta gaiiHi which, for certain param- 
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FIG. 2. (Color online) Density profiles n(z) of the r' = 0 soliton for (a) ^ = 0 and (b) Q — 'k as a function of Sdd- The band 
(grey) of instability is bounded by the onset of the phonon and roton instabilities (dashed and solid lines, respectively). Inset: 
density prohle for no dipoles (sdd = 0) and close to the instabilities (sdd = 0.9,1.45). Note that the unit of length ^ is itself a 
function of edd and 0. 


eter regimes (red regions in Fig. [^b)) can touch zero 
energy, signalling the unstable growth of finite momen¬ 
tum modes, i.e. the roton instability. When deep in the 
ID regime (a <C 1, as employed here) the roton instabil¬ 
ity arises for < 0 [48] [Fig. Bb)] ; however, for a ^ 1 
the roton instability shifts to ^ > 0, as predicted in Ref. 

m- 


IV. DARK SOLITON SOLUTIONS 

Having established the stability of the homogeneous 
quasi-ID dipolar BEG, we move on to the dark soliton 
solutions themselves. It is known that the purely local 
GPE is integrable and supports a family of dark soliton 
solutions for repulsive contact interactions (^ > 0) [50l 

m], 


V’s(^d) = 


■ , p(z-Z) .v' 

/5tanh -— -h i- 




(4) 


Here [3 = y/l — and Z{t) = -\- vt^ where zq is 

the soliton’s initial position and v its speed. Stationary 
solitons have a node of zero density and a phase slip of tt, 
whereas a = c soliton has no phase or density contrast 
from the background fluid. The soliton energy decreases 
with increasing speed m, leading to the analog of a 
particle with negative effective mass [52] . Meanwhile the 
density minimum of the soliton scales as Umin/^o = 

0 . 

Noting that dark solitons are stationary solutions in a 
moving frame, we numerically seek solutions to the GPE 
§ with a Galilean boost term ihvdzip in the Hamil¬ 
tonian. Discretizing t/; on a ID spatial grid (spacing 
dz = 0.1^) we minimise the discretized Hamiltonian with 
respect to changes in -0 using a bi-conjugate method [^ . 
Starting with the non-dipolar dark soliton solution Q, 
this method leads to convergence to the required dipolar 
soliton solution, within the entire speed range 0 < v < c. 


The numerical box (up to ±800^) is large enough to 
mimic the infinite limit. 

The numerically-obtained dark solutions propagate 
with constant speed and permanent form (no dispersion 
or radiative losses) when simulated within the lab-frame 
GPE [55], confirming their soiitonic character. Dark soii- 
tons require the mean-fieid potentiai of the homogeneous 
system to be net repuisive [561157], i.e. /io > 0, the same 
as the condition for phonon stabiiity. 

We iiiustrate the soiutions through the v = 0 soiiton: 
Eig. depicts its spatiai density profiie n{z) as a func¬ 
tion of Sdd foi* the limiting angles (a) 6> = 0 and (b) 
0 = it/2. The former is indicative of the general be¬ 
haviour for 0 < Om, and the latter for 0 > O^n (for 0 = 0^ 
one recovers the non-dipolar soliton throughout). The 
soliton size is characterised by the dipolar healing length 
^ (hence motivating this choice of units). Away from the 
unstable band of the profile approximates the tanh- 
squared density profile of the non-dipolar soliton. As the 
phonon instability is approached (for ^ > 0) the profile 
diverges in width; this is related to a cancellation be¬ 
tween the local interactions arising from the explicit van 
der Waals interactions and an implicit local contribution 
to the dipole-dipole interactions [59]. Meanwhile, as the 
roton instability is approached, prominent density ripples 
form about the soliton core (see inset). These finite-/c 
corrugations arise due to the mixing of the roton with 
the soliton state, with analogous ripples arising for vor¬ 
tices [59f[6T] . With increasing speed, the soliton depth 
decreases in line with the non-dipolar soliton, and the 
above qualitative behaviour remains, albeit with weaker 
ripples as the speed is increased. The soliton phase pro¬ 
file is insensitive to e^d and 0. 

The soliton modifies the mean-field dipolar potential 
^id(^) in a non-local manner. We can glean some in¬ 
sight into this by considering the behaviour at long range. 
Expanding the dipolar pseudo-potential Eq. around 
infinity gives Uib{u) = Uo{4:\u\~^ — 24\u\~^ + ...) [62] . 
Then, taking the non-dipolar soliton as an ansatz, the 
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asymptotic form of 4 >id follows as, 





— z\^J ’ 


( 5 ) 


where zq is a short-range cut-off to account for the 
asymptotic behaviour of the mean-field dipolar poten¬ 
tial. The two terms represent the background and soli- 
ton contributions to 4 >id, respectively. The 1/z^ decay 
of 4 >id demonstrates that the soliton appears as a lo¬ 
calized giant dipole, a key result of our work. This is 
intimately connected to the exponentially-fast decay of 
the soliton density profile; at large distances the soliton 
profile scales as n{z) ^ (no — In contrast, a 

slow power-law-decaying density profile would render the 
object as an extended dipole at all scales, as for vortices 
in 2D [59] . The negative sign on the soliton contribution 
suggests that the soliton may be viewed as a collection 
of anti-dipoles superposed on a homogeneous background 
of conventional dipoles. This observation will help us be¬ 
low to interpret the modified form of the soliton-soliton 
collisions as due to the interaction of these anti-dipoles. 


V. SOLITON-SOLITON COLLISIONS 

We can expect the non-locality of the solitons to mod¬ 
ify their interactions, which we examine next through 
their collisions. For simplicity we take the collisions to 
be symmetric, i.e. same incoming speed for both soli- 
tons. The corresponding analytic two-soliton solution to 
the non-dipolar ID GPE is [5], 

= {F{z,t)/G{z,t)}e-x.p{-iiJLt/n) ( 6 ) 


with, 

F{z,t) = 2^/no [(1 — 2v^/c^)cosh{2v 131/^) 

— (v/c) cos\i{2l3z/£^) + i2{v/c)l3sin\i{2vl3t/£^)] , 

G{z,t) = 2cosh{2vI3t/^) F 2{v/c)cosh{2^z/^). 

Conventionally, solitons are known to repel during their 
collision with two regimes of dynamics: for low incoming 
speeds v < 0.5c the solitons reflect at short-range, with 
a sudden large repulsion akin to hard-sphere collisions, 
while for higher speeds v > 0.5c the solitons transmit 
dZlEs]. We focus on slow collisions v = 0.1c, comment¬ 
ing on the wider speed dependence below. We adopt 
fixed dipole parameters corresponding to ^^^Er [35], i.e. 
^dd = 0-4. Eigure 3 shows some example soliton colli¬ 
sions, simulated via the ID dipolar GPE [55], for differ¬ 
ent values of polarization angle 0. The solitons emerge 
from the collisions with unchanged profile and no radia¬ 
tive loss, further supporting their solitonic character. At 
(b) the magic angle Om, one recovers the non-dipolar dy¬ 
namics in which the solitons “bounce” at short-range [63] . 
Eor (a) 0 < the soliton repulsion is significantly en¬ 
hanced, turning them at increased separation, indicative 
of the role of the non-local interactions. Here the dipoles 



FIG. 3. (Color online) Collisions of dipolar dark solitons with 
incoming speeds v = 0.1c for Cdd = 0.4 and polarization angles 
(a) ^ = 0, (b) ^ = Om and (c) 0 = 7r/2. 


are net attractive but the solitons themselves (which, re¬ 
call, behave as anti-dipoles) are effectively net repulsive. 
Vice versa, for (c) 0 > 0^ the dipole-dipole interactions 
introduce an attraction between solitons; however, the ef¬ 
fect of this is barely visible due to the dominance of the 
hard-sphere-like repulsion at short-range. 

The non-local repulsive/attractive modifications to the 
soliton collisions become more pronounced with increas¬ 
ing Sdd and a due to the increasing strength of the ef¬ 
fective ID dipole-dipole interactions. Meanwhile, they 
become less pronounced with increasing soliton speeds; 
this is due to the reduced mean-field potential generated 
by faster and shallower solitons, as evident from the (3 de¬ 
pendence in Eq. <§• While the elastic collisions in Eig. 
3 are representative of most of the stable (6>, 5dd)-space, 
collisions become inelastic close to the phonon instability 
or roton instability; then sound waves are shed during the 
collision, causing a slight (few percent) increase in their 
outgoing speeds. 


VI. TWO-SOLITON INTERACTION 
POTENTIAL AND BOUND STATES 


We can further understand the soliton-soliton interac¬ 
tion through an effective particle-like interaction poten¬ 
tial. We calculate the soliton’s energy, relative to the 
background, as. 


-Etot — -Eq + Edd-i 


( 7 ) 


where. 


Ed 



'F 

2m 


\d.F + l(\F 



d^. 


(8) 


is the non-dipolar energy (sum of kinetic and van der 
Waals interaction energies) and. 


^dd 


/ 


1 

2 


^id\F 


( 9 ) 


is the dipolar contribution. Eor two solitons at zi and 2:2 
(defined by their density minima), with separation q = 
\zi — Z 2 \, we define the interaction energy as. 


V{q) = Etot{zi,Z 2 ) - Etot{zi) - Etot{z 2 )- (10) 
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FIG. 4. (Color online) (a) Soliton-soliton interaction poten¬ 
tial y as a function of separation q for Sdd = 0.4 and Q — 'k 
(estimated using the non-dipolar two-soliton solution in the 
low speed limit with arbitrary value v — 0.001c). The local 
minimum indicates a bound state. The non-dipolar Vb and 
dipolar Vdd contributions are indicated, (b)-(d) Correspond¬ 
ing CPE dynamics starting from two stationary solitons a 
distance go apart. For (a) go = 4^, the solitons repel, while 
for (b) go = 5^ and (c) go = 10^, bound state oscillations are 
evident. 


We estimate T^(g) semi-analytically based on the non- 
dipolar two-soliton solution with small incoming speed 
(for larger speeds, the effective masses of the solitons 
change considerably during the collision, complicating 
this particle-like picture). Figure |^a) shows the soliton- 
soliton interaction potential, again for ^^^Er parameters. 
The non-dipolar contribution Vb (dashed red line) dom¬ 
inates only at short-range, consistent with the repulsive 
bouncing of two non-dipolar solitons [Fig. 3(b)]. For 
^ = 0, the dipolar interaction potential Vdd (black dot- 
dashed line) is repulsive and non-local, consistent with 
the bouncing at increasing separation observed in Fig. 
3(a). However, for d = 7r/2, Vdd (blue dotted line) is 
attractive at long-range. The playoff between this non¬ 
local attraction and the short-range repulsion from Vb, 
conspires to form a total potential (solid black line) anal¬ 
ogous to the Lennard-Jones inter-atomic potential: an 
energy minimum at finite g, with rapidly increasing po¬ 
tential at short-range and slowly increasing potental at 
larger range. This raises the prospect of forming a two- 
soliton bound state, analogous to a diatomic molecule. 

The possibility of supporting a two-soliton bound state 
is probed through CPE simulations starting with two 
stationary solitons a distance go apart, shown in Fig. 
[^b)-(d). For sufficiently small go (a), the solitons ini¬ 
tially feel a strongly repulsive short-range interaction, 
repel and acquire sufficient kinetic energy to escape to 
infinity. For greater separations, e.g. (b) and (c), the 
solitons initially have negative potential energy and are 




FIG. 5. (Color online) Collision of two counter-propagating 
bound states, with speeds (a) v — 0.01c and (b) v — 0.5c 
In both cases, the bound states emerge unscathed from the 
collision, with the only net effect being a phase shift of the 
outgoing waves (for the high speed case this phase shift is 
only just visible). 


restricted to execute oscillations in g about the potential 
minimum, analogous to vibrational modes of a diatomic 
molecule. For large initial separations go ^ ^min, tbe 
solitons experience a weak gradient in V (g) and undergo 
slow oscillations. 

Well away from the roton instability and phonon insta¬ 
bility, these bound-state oscillations persist ad infinitum 
(within the zero-temperature CPE). However, close to 
these instabilities, repeated sound emission during the 
collisions leads to the counter-intuitive situation where 
the solitons increase their oscillation amplitude and ul¬ 
timately escape the bound state by losing energy; this is 
related to the negative effective mass of the soliton and is 
analogous to anti-damping effects in external potentials 

m- 

We further investigate the properties of the bound 
state by considering its interaction with another bound 
state. Figure depicts the collisions of two counter- 
propagating bound states at (a) low and (b) high speed. 
These moving bound states are formed by repeating the 
above method for forming bound states, but where both 
initial solitons have the same non-zero speed. First, it is 
worth observing from the plot that the bound states are 
stable to centre-of-mass motion at constant speed. For 
both low and high speed, the two bound states emerge 
unscathed from the collision, with their original speeds 
and with no radiative losses. The only net difference is 
that the outgoing bound states both feature a transla¬ 
tional offset, termed a phase shift. This is considerably 
larger for the low speed case, and is just visible in the 
high speed collision. At low speed, the bound states ap¬ 
pear to bounce off each other, just like in the collision of 
two slow dark solitons [Fig. 3], while at high speed the 
bound states appear to pass through each other, again 
analogous to the corresponding behaviour of two fast col¬ 
liding dark solitons [5]. Note that for the high speed case, 
the solitons in each bound state appear to move in paral¬ 
lel; this is simply because the period of the bound state 
oscillation is considerably longer than the timescale of 
the figure and the collision. We have repeated these sim- 
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ulations over a wider rage of incoming speeds, and find 
the same qualitative soliton-like behaviour throughout. 
These results demonstrate the striking property that the 
bound states themselves behave like solitons in their in¬ 
teractions with other bound states. 

Furthermore, analysing the collision between a bound 
state and a single incident dark soliton shows the same 
behaviour, with the single soliton and bound state emerg¬ 
ing unscathed from the collision, barring a phase shift. 


VII. CONCLUSIONS 

In summary, we have studied self-trapped non-local 
dark solitons supported within quasi-lD dipolar BECs. 
The solitons acquire modified profiles, including ripples 
associated with roton excitations. The solitons approxi¬ 
mate giant localized dipoles, and have non-local soliton- 
soliton interactions, controllable through the direction of 
polarization of the dipoles. When attractive, and in com¬ 
bination with the conventional short-range repulsive in¬ 
teraction, unconventional dark soliton bound states can 
be realized. These bound states are stable to centre-of- 
mass propagation at constant speed. Remarkably, they 
act themselves like solitons during collisions, emerging 
with unchanged form and speed. 

Some of these effects are analogous to predictions 
for vortices in two-dimensional dipolar condensates m- 
There, vortices bear similar ripples about the core, with 
these ripples being a common manifestation in dipolar 
condensates in the vicinity of a sharply-varying density 
profile. However, the different density profiles of solitons 
and vortices lead to signficantly different results. For a 
vortex, the profile scales asymptotically as 1/r^ back to 


the background density; this is sufficiently slowly vary¬ 
ing that it renders the vortex as an extended dipole at 
all scales. In contrast, the asymptotic part of the soli- 
ton density profile decays exponentially fast to the back¬ 
ground, such that the soliton resembles a localized dipole 
at long range. These differences, in turn, lead to differ¬ 
ent functional forms of their interaction potentials: for 
solitons the interaction potential scales with the dipo¬ 
lar form 1/r^ at long range, while for vortices the most 
slowly decaying terms go as 1/r^ and ln(r)/r^. For vor¬ 
tices this is in addition to the usual long-range, isotropic 
interaction arising from hydrodynamic effects. 

The above soliton behaviour occurs robustly across a 
wide parameter space, accessible to dipolar BEC exper¬ 
iments |3Tll33l [ 35 ] and with current dark soliton gen¬ 
eration and imaging capabilities [iSHn]- These dipolar 
dark solitons extend the range of physics of dark solitons 
as mesoscopic probes of quantum physics m to include 
the interplay with magnetism. Moreover, the effectively 
instantaneous non-local soliton interaction, and its ca¬ 
pacity to be experimentally tuned at will, offers intrigu¬ 
ing possibilities for the controlled study of non-locality 
in complex networks [25], soliton gases [66| and super- 
solitons 

Note Added: After completing this research we became 
aware of the manuscript [68| which obtains results related 
to ours. 
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